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Abstract 



We study the quantum tunnelling of a very complex object of which only 
part is coupled to an external potential ( the potential barrier ). We treat 
this problem as the tunnelling of a particle (part of the system affected by the 
potential) coupled to a "heat bath" ( all of the remaining degrees of freedom 
) from the point of view of the wavefunction of the whole system. We obtain 
an effective scattering equation for this wavefunction which resembles the 
integral equation for the wavefunction of the particle in the absence of the 
bath. We show that the effect of the interaction is to renormalize the mass 
and the potential of the problem and that these new parameters can be solely 
described in terms of the spectral function of the excitations of the system. 

PACS numbers: 03.65.-w, 05.30.-d, 05.40.+j, 03.65.Nk 
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The problem of the quantum mechanical tunnelling of non-isolated particles has been at- 
tracting a great deal of attention for over one decade now. People have been maily interested 
in the so-called dissipative quantum tunnelling which turns out to be of crucial importance 
to problems involving either macroscopic quantum tunnelling or macroscopic quantum 
coherence .Although the mathematical details can vary from one problem to the other, 
the basic idea behind all of them is always the same. One tries to mimic the influence of 
a very complicated set of degrees of freedom by a set of non-interacting harmonic modes 
with an appropriate spectral function previously chosen . The next step is to trace out 
these harmonic modes from the dynamical or equilibrium density operator of the whole sys- 
tem, depending on the problem one wants to attack, and use the so-called reduced density 
operator of the system in order to obtain the desired information. 

The problem we want to face in this work is not of the same kind as the above-mentioned 
ones. What we are now interested in is the quantum tunnelling of a very complex system, 
say, a huge molecule or a very heavy nucleus, of which only a small set of variables can 
couple to an external potential. We prepare our complex system as a wave packet moving 
towards that barrier and ask the following question: "How do the extra modes ( the ones not 
feeling the potential barrier ) affect the transmission coefficient of the complex system as a 
whole?" 

As an example, suppose we want to study the tunnelling of a very complex nucleus 
through a potential barrier of purely electromagnetic origin. In this case only the protons 
will feel the presence of the potential barrier although they will still be coupled to the 
neutrons via strong interaction. If there was no external potential, the total momentum of 
the system would be conserved. However, the presence of this barrier breaks the translation 
invariance of the problem and the nucleus will be either reflected or transmitted. This takes 
place through the motion of the protons which on their turn influence the motion of the 
neutrons. Therefore, the latter will ultimately affect the motion of the center of mass of the 
system and we can test it by studying the effective reflection and transmission coefficients 
of the nucleus whole. 

In order to simplify matters and to keep away from intricate microscopic models we 
shall model our complex system as a particle (the relevant degrees of freedom) coupled to a 
bath ( the complex set of internal excitations of the system) subject to a potential barrier 
acting only on the particle "coordinate" .The bath on its turn will be described by a set of 
non-interacting oscillators and due to our hypothesis of its high complexity we shall also 
define its spectral function and "internal temperature" . 
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Therefore, from now on, we shall be employing the same terminology as for the general 
problem of quantum dissipation bearing in mind that it should be taken seriously only in a 
very strict sense. We shall emphasize this remark whenever it is necessary. 

Although stated in this way this problem resembles a lot the one of the transmission of 
a dissipative particle they are not quite the same. For the latter the boundary condition 
must be entirely different from the one we shall be using in the present work. Here, the 
"system-plus- reservoir" move together in the same state of motion towards the barrier. It 
is clear that the presence of the external potential may cause "internal dissipation" due to 
the flow of energy from the particle to the bath. Nevertheless, the total energy of the system 
will be conserved. 

It will be shown that the system tunnels as a particle with renormalized mass through a 
temperature dependent effective potential barrier. Moreover, the parameters of this effective 
scattering theory will be dependent only on the form of the spectral function of the internal 
excitations of the system. 

Our approach in this article is quite different from the usual one for dealing with the 
quantum dynamics of non-isolated systems. Instead of treating the time dependent problem 
using path integrals we will use a time independent approach which resembles the problem 
of scattering by an external potential, in other words, we will write the time independent 
Schrodinger equation for + 1 bodies in the approximation where the internal degrees of 
freedom act like a heat bath for the particle that feels the external potential. 

Our starting point will be the well-known hamiltonian of a particle coupled linearly to 

harmonic oscillators which is described as follows [|],^, 

where p and x, pi and Xi, are the momentum and the coordinate of the particle and the 
oscillators, respectively, m is the "bare" mass of the particle, Ui and rrii are the frequency and 
the mass of the i^^ oscillator, respectively, Ci is the coupling constant between the particle 
and the oscillators and V{x) is the external potential applied to the particle. 

It was shown by Hakim and Ambegaokar [^] that the hamiltonian (1) can be diagonalized 
by a unitary transformation when there is no external potential {V = 0). In this case we 
obtain a hamiltonian where the center of mass of the whole system, particle plus oscillators, 
is decoupled from a new set of noninteracting oscillators. The diagonalized hamiltonian can 
be written as [01, 
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H' = U^HU = + y hVL.a\a. (2) 

2M j 

where f/ is a unitary transformation, P is the momentum of the center of mass of the system 
and Af is the renormahzed mass of the system which is given by, 

N 

M = m + J2N (3) 



with 



/^i = — —■ 

rrijUJj 

The operators Oj and aj are the creation and annihilation operators for the new harmonic 
modes with frequency Qj which are solutions of the equation 

y ^ = 1 

The connection with the dissipative problem can be made if we introduce the spectral 
function of the oscillators, 

Ji^) = E ^Siu; - ^0 (4) 
which in the ohmic regime is described by ||^, 

J{uj) = rjuQ^uJc — uj) 

where r] is the damping coefficient and uJc is the cut-off frequency of the bath. However, one 
should notice that this kind of assumption for the behavior of the spectral function is not 
appropriate for the present problem. 

Observe that in infinite systems the total mass M, eq. (3), may diverge depending on 
the behavior of (4) at low frequencies. This shall not be the case for the sort of systems in 
which we are interested in this paper. We shall restrict ourselves to finite systems and these 
are bound to present a gap in their spectrum of excitations. In other words, the spectral 
function (4) will always have a natural low frequency cutoff. 

In order to proceed further we include the external potential in the problem. Following 
the route established in ref. |0 we include a new term in (2) which is given by U'^VU where 

f/ = e-^S^''^^"^''^^ (5) 
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with 



N I, n^/'^ ( N ,, , ,2 



1=1 



1/2 



r \.2 (6) 



or in terms of the coordinates of the new oscillators we have 

U'^VU^vi^-Y.^ ^2hf,jnj , (7) 

where X is the position of the center of mass of the system. 

Equations (6) and (7) show that the center of mass is coupled to the renormalized 
oscillators via the external potential. Suppose now that the potential has a finite range, 
in other words, it is a potential barrier. Moreover, since there are no modes of very low 
frequencies in the problem, there will always be a region for X sufficiently far from the 
origin where (7) vanishes. Therefore, outside the scattering region we have a free particle 
with mass M and momentum hq and a set of decoupled harmonic oscillators in states 
described by occupation numbers {Nj}. This is a scenario for scattering since we already 
know the asymptotic states of the problem. This is not revealed in our original hamiltonian 

(!)• 

The most natural way to solve this problem is to rewrite the time independent 
Schrodinger equation in the representation of the momentum of the particle and the Fock 
space of the oscillators. Adding (7) to (2) and writing the resulting equation in this repre- 
sentation together with the above-mentioned incoming condition, we easily get 

^E{k,{nj})^S{k-q)ll6n„N,- E W7r2 2^ t,o ( 

where is the wave function oi N+1 bodies with total energy E — + Z^j hQjNj, V{k) 
is the Fourier transform of the potential and finally, 

G{k - ^'){n,},K} M\U\k)U{k')\{n[} > (9) 

is the thermal part of the kernel of equation (8) where U is defined in (5). Observe that the 
effect of the external potential is completely separated from the coupling to the oscillators 
which is present only in (9). This remarkable fact will allow us to study the problem 
independently of the form of the external potential. 

Notice that when the interaction in turned off, that is, Cj = 0, we have U = I where / is 
the identity operator. Therefore, G{k — k')^ni},{n'.} — ^{ni},{n'.} ^-nd the solution is given by 
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'^E{K{ni})^^l^^{k)5{r^,UN,} (10) 

where obeys the "bare" scattering equation 

(fc) = 5{q -k)- I E ^(^ - (11) 

■^{1^ -q ) k' 

as expected. 

The solution of eq. (8) can be very difficuh depending on the shape of the external 
potential (as in the problem without dissipation, eq. (11)). The presence of the kernel (9) 
creates new problems since it contains all induced transitions between the states of the 
oscillators due to the recoil of the core particle. We would like to mention that there is a 
formal solution for (8) for the case of local interactions {V{x) — aS{x)) but the solution is 
not illuminating due to its complexity 

Our approach here is to suppose that there is a huge number of normal modes and they 
are in thermal equilibrium. Notice that we are by no means planning to model the internal 
modes of our complex system as an infinite reservoir. By thermal equilibrium we only mean 
a fixed average energy which allows us to define a parameter, the "internal temperature" , 
we shall loosely refer to as simply temperature. What we aim at with this approximation 
is to have the system of oscillators being only very weakly influenced by the recoil of the 
particle. Although this hypothesis is more reasonable for high temperatures we shall assume 
it is also the case at very low temperatures. We then replace (9) by its average value in 
thermal equilibrium which is given by, 

G{k){n,},{nj} = hmUnj} < G{k) >t (12) 

where 

tvR (e-^"'^G(k)) 
< G{k) >T= ^ , \ (13) 

where tr^ means the trace over the reservoir coordinates, Hji is the bath hamiltonian given 
by the last term on the r.h.s. of (2) and /3 = {KbT')~^ where T is temperature {Kb is the 
Boltzmann constant). 

If we substitute (12) in (8) we find the effective scattering equation for the center of 
mass, 

U^) = b{k - g) - ^ J E ^'^^(^ - ^) U^) (14) 

-mv^ ~ q ) k' 
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where 



V'f^{k, T) - V{k) < G{k) >T (15) 

is the effective potential felt by a "dressed" particle with mass M . Now the kernel of the 
scattering equation depends on the temperature and the coupling. The solution of (14) will 
provide us with the reflection and transmission coefficients for the scattering of the particle 
in the presence of a large set of internal degrees of freedom. It is worth mentioning that the 
wave function ip is not the wave function for the bare particle but it is the effective wave 
function for the particle dragging the surroundings in its motion. 

We can evaluate the thermal part of the kernel given in (13). We can show, after some 
straightforward, but lengthy, algebra that 

< G{k) >r= e-r(^)*' (16) 

where 

r(T) = Vir|(l + 2n,), (17) 



and 



is the Bose occupation number. 

If we substitute (16) in (15) and transform the problem to real space we find that the 
effective potential is written as 

V^^^{X, T) = J dX' V{X') A(X - X\ T) (18) 

where 

and therefore if the interaction is turned off (F 0) we obtain the result that the effective 
potential is exactly the bare potential {Limr^oA{X,T) = 6{x)). 
First, we note that the area under the potential is conserved, 

J dxV^f^ix) = J dxV{x) (20) 



7 



as we can easily check. Moreover, it is easy to show that the height of the potential is a 
monotonic decreasing function of r(T'). Since the area is conserved we conclude that the 
width of the potential must increase. 

As an illustration of our results consider the simple case of a gaussian potential, 

l/(a;) = yoexp|-^|. (21) 

It is very easy to obtain from (18) that the effective potential is given by. 

Observe that the effect of interaction was to decrease the height of the potential and increase 
its width accordingly (see equation (20)). Indeed we have. 



V^ff{0) a \ a?' 

The effect on the transmission coefficient is readily obtained. Prom energy conservation 
we have. 



nq= ^J2M{E-Et)\ 



hK = ^2M(K//(0)-£; + £;t), (24) 

where E is the total energy of the system whereas Et = Y^j hQjNj stands for its internal 
energy. For temperatures such that na » 1, we can show that the transmission is given 

by, 

|t|2 ^ 16 ( e-^-<->eff_ (25) 
+ J 

The exponent of the transmission coefficient in (25) can be easily analyzed if we write it 

as 



Ak<x >eff^ 4^2MVoa^^Vl + 4.r^ - (1 + 4r2) / (26) 

where r = T/a , f = Ex/Vq and E^ = E — Et which is the energy of the center of mass of 
the system. Observe that we have two effects to take into account in (26). The increase of 
the effective mass, M, of the particle and its dependence on the parameter r. 

The mass dependence has an obvious influence on the transmission coefficient because 
(26) monotonically increases with it. Therefore the transmission will steadly decrease as a 
function of the effective mass which is an expected result. 
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On the other hand, the dependence of the transmission coefficient on the combination of 
the effective height and width of the potential is not so simple. For a given spectral function 
one can compute the minimum effective range of the potential by using (17) for T = O.This 
provides us with the minimum value of r since F is an increasing function of temperature. 
Now, it is worth noticing that (26) has a maximum at 1/4/ and vanishes at 1/2/ when 
/ << 1 . Therefore one can say that as the coupling is switched on there is a reduction of 
Itp as compared to its non-interacting value. Furthermore, as one increases the temperature 
there is first a reduction of the transmission coefficient before it starts to increase to its 
maximum value. Actually we are assuming that in the limit of small / expression (25) is a 
good approximation for any value of r^m < r < 1/2/. Actually, this is true because in this 
range na > 1. 

One can therefore conclude that the effect of the coupling to the internal degrees of 
freedom of the system tends to reduce its rate of transmission through potentials barriers. 
Actually there is a competition between this effect and the one related to the decrease of the 
potential barrier seen by the system as its center of mass kinetic energy is increased. The 
latter wins when 1/4/ < r < 1/2/ and this causes the transmission coefficient to increase 
in this region of energy. 

This result is in agreement with those obtained in the study of dissipative quantum 
tunnelling and coherence |l|,0] which say that the coupling of the form (1) to a heat bath 
inhibts quantum mechanical effects. Actually at this point one should be really tempted 
to draw some conclusions about dissipative quantum transmission from the results we have 
got in this paper. However, although we think that the final results for the latter will 
not differ much from the present ones, the reader should be warned that the boundary 
conditions appropriate to that problem are completely different from the one we used here 
and, in fact, it is not even clear to us whether model (1) in its present form should be the 
one to describe this phenomenon. The application of our present results with a spectral 
function J {u) = acu" should be very carefully interpreted since we are treating only finite 
environments in this work. It would only make sense for very large n since in this way we 
should be trying to simulate a gap in the spectrum of excitations of the system. 

In conclusion, we have derived an effective scattering equation for a particle coupled 
to a large, but finite, set of oscillators (internal degrees of freedom) which resembles the 
scattering equation for an isolated particle. We found that the main effects in the scattering 
problem are the increase of the mass due to the dragging of the oscillators and a change in 
the scattering potential which depends on the temperature and J (tu) (through the coupling 
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constants of the problem). We showed that the results for high and thin barriers can be 
obtained independently of the specific form of the bare external potential and they imply 
in the reduction of the transmission coefficient for tunnelling at low center of mass energies. 
The main assumption in our derivation is that the oscillators retain their thermal equilibrium 
during the scattering process. 

A. H. Castro Neto gratefully acknowledges M. P. Gelfand and A. J. Leggett for useful 
comments, F. Guinea for illuminating discussions, Conselho Nacional de Desenvolvimento 
Cientifico e Tecnologico, CNPq (Brazil), for a scholarship. A. O. Caldeira also wishes to 
acknowledge the support from the Conselho Nacional de Desenvolvimento Cientifico e Tec- 
nologico, CNPq (Brazil). This work was partially supported by the National Science Foun- 
dation through grants no. PHY 89-04035 and DMR 92-14236. 



10 



REFERENCES 
[1] A.O.Caldeira and A.J.Leggett, Ann.Phys. (NY) 149, 374, (1983). 

[2] A.J.Leggett, S.Chakravarty, A.T.Dorsey, M. P. A. Fisher, A.Garg, and W.Zwerger, 
Rev.Mod.Phys. 59, 1, (1987). 

[3] A.O.Caldeira and A.J.Leggett, Physica (Amsterdam) 121 A, 587, (1983); 130 A, 374(E), 
(1985). 

[4] R.P.Feynman and F.L.Vernon, Ann.Phys. (NY), 24, 118, (1963). 
[5] R.Bruinsma and Per Bak, Phys. Rev. Lett, 56, 420, (1986). 
[6] V.Hakim and V.Ambegaokar, Phys.Rev.A 32, 423, (1985). 
[7] A.H.Castro Neto, unpubhshed. 



11 



